function [mu_U, L_U] = rnd_user_hyperparams(V)
[D,N] = size(V);

beta0 = 1;
nu0 = D;
mu0 = zeros(D,1);
W0 = eye(D);

V_ = mean(V,2);
S_ = fcov(V,2);

beta0_ = beta0 + N;
nu0_ = nu0 + N;
mu0_ = (beta0*mu0 + N*V_)/(beta0+N);
W0_ = inv(inv(W0) + N*S_ + (beta0*N/(beta0+N))*(mu0-V_)*(mu0-V_)');

W0_ = (W0_+W0_')/2;
L_U = wishrnd(W0_, nu0_);
Li_U = inv(beta0_*L_U);
mu_U = rnd_mvn(mu0_, Li_U);
